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The Korringa ratio, /C, obtained by taking an appropriate combination of the Knight shift and 
nuclear spin-lattice relaxation time, is calculated at finite temperature, T, in the three-dimensional 
electron gas model, including the electron-electron interaction, U, and non-magnetic impurity scat- 
terings. /C varies in a simple way with respect to U and T; it decreases as U is increased but increases 
as T is raised. However, K. varies in a slightly more complicated way with respect to the impurity 
O i. scatterings; as the scattering rate is increased, tC increases for small U and low T, but decreases for 

large U or high T regime. This calls for a more careful analysis when one attempts to estimate the 
^0 ' Stoner factor from fC. 
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£ ■ I. INTRODUCTION 

c/3 ; 

The Korringa ratio, /C, given by K.= 1/(T\TK 2 ), where T\ is the nuclear spin-lattice relaxation time, K the Knight 

S shift, and T is the temperature, is a useful concept [n|, in that when K. of a material is larger (smaller) than unity, it 

is interpreted as a strong indication that the material is antiferromagnetically (ferromagnetically) correlated. Moriya, 
i-rt ! some forty years ago, showed that /C is decreased when there is electron-electron interaction, U, in a three-dimensional 
(3D) free electron gas ||. Underlying physics is that K 2 increases more rapidly than \jT\T for a ferromagnetically 
correlated system as U is increased. It was later extended to include the non-magnetic impurity scattering effects 
by Shastry and Abrahams ||. They showed that disorder enhances the Korringa ratio so that the impurity effects 
should be included in addition to the electron-electron interaction when one analyzes the ratio, eg, when one wants to 
t-H , estimate the Stoner factor of a material from /C. These calculations were done at zero temperature. Since experiments 
are performed at finite temperature, we wish to extend their calcuations to the finite temperature, and understand 
how K. is changed as the temperature, T, the electron-electron interaction, U, and the impurity scattering rate, 1/r, 
JLJ ■ are varied. 

K and 1/Ti, which determine K,, can be expressed in terms of the spin susceptibility, x(q, uj), of the conduction 
electrons. Fulde and Luther calculated the impurity effects on the spin susceptibility of almost ferromagnetic metals 
f^i ' by including isotropic non-magnetic impurity scatterings and Dirac ^-function electron-electron interaction at zero 
temperature, in the small frequency and the small wavevector limit B. Their result was utilized by Shastry and 
Abrahams to calculate K. of almost ferromagnetic dirty metals. In the present work, we reformulate the spin suscep- 
tibility with the Matsubara Green's function to consider the finite temperature effects on the Korringa ratio at the 
\~ | presence of the electron-electron interaction and impurity scatterings. The calculated K, T%, and /C agree in most 

cases with those of Shastry and Abrahams: both K and \jT\T increase when U or 1/t is increased, and K, decreases 
as U is increased. One interesting observation, different from their results, is that /C increases as 1/t is increased only 
when both U and T are small. When U is comparable with U cr , or T is with ep , /C decreases as 1/r is raised. U cr 
Q , is the critical value of the electron-electron interaction at which the ferromagnetic instability occurs, and ep is the 
!LJ Fermi energy. Consequently, K, vs. U curves for several different 1/r values, at a fixed T, cross each other. These 
. ,— i \ crosses also occur in K. vs. T at a given U among different values of 1/r. It has been believed, since the Shastry 
and Abrahams work, that the impurity scatterings enhance /C while the electron-electron interaction suppresses it. 
The present result that impurity scatterings can enhance or suppress depending on U or T requires that the above 
interpretation should be modified. This calls for a more careful analysis of experimental data when one wishes to 
estimate the strength of ferromagnetic correlation of a material. 

The remainder of the present paper is organized as follows: In Sec. II, we present the Matsubara formalism to 
include the electron-electron interaction and impurity scatterings in calculating x(q> w ) a t hnite temperature. This 
was carried out in the Matsubara frequency and then analytically continued to the real frequency. In Sec. Ill, K, 
\fT\T ', and /C are calculated by numerically integrating the obtained equations, as U, 1/r, and T are varied. These 
results are compared with the free electron gas model as a reference, and with previously reported calculations. 
Finally, Sec. IV is for the summary and concluding remarks. 



II. FORMALISM 



The Korringa ratio is obtained from the Knight shift and nuclear spin-lattice relaxation time as JC— l/(TiTK 2 ), 
and K and 1/Ti are given in terms of the spin susceptibility as 
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Calculating the Korringa ratio, therefore, amounts to calculation of x(q> w ) including the electron-electron interaction 
and the impurity scatterings at finite temperature. This can most conveniently be done in terms of the Matsubara 
Green's function in the imaginary frequency to obtain x(q,iui). It is then analytically continued to yield x(q, w) in 
the real frequency by substituting iu — > u + iS, where 6 is a positive infinitesimal. We, therefore, reformulate the 
Fulde and Luther result with the Matsubara method for spin susceptibility and construct K and (TiT) -1 . The spin 
susceptibility including the vertex corrections from the Coulomb interaction and electron-impurity scatterings can be 
expressed in terms of the Feynman diagram as in Fig. 1(a). Here, the solid lines denote the renormalized electron 
Green's function G(k, ip), and T denotes the vertex function, where ip and iu are, respectively, the fermion and 
boson Matsubara frequencies, and k and q are the correspinding wavevectors. T includes both the impurity and the 
Coulomb interaction as shown in Fig. 1(b), but G does not include the renormalization from the electron-electron 
interaction as in the previous works j|||. The neglect of U in calculating G may be justified because Schrieffer and 
Berk showed that this model reproduces the long wave length static susceptibility which would be obtained from a 
more complete calculation including the renormalizations due to U ||. 
r, represented in Fig. 1(b), can be written as 
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Here, (3 = l/fc^T, TV the single-spin density of states (DOS) at the Fermi level, and G(k, ip) = l/[ip+ ^sgn(p) — £k], 
where sgn(p) is 1 if p is positive or zero and is —1 otherwise, and £k is energy measured from the Fermi surface, or 
£k = fi, 2 k 2 /2m — €f- X can be written as in Ref. B in Matsubara form as 

2 
X(q, iu) = -n% -^^] G(k, ip) G(k + q, ip + iu) T(q, ip, iu) (3) 
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sgn(p)y/ip+/j, + sgn(p+u)^i(p + u)+fi + 
sgn(p)^/ip+jl + sgn(p+u) y/i(fi+u)+p, - A 



Here, we put ks = ti = 1, /is the Bohr magneton, m the electron mass, e q = q 2 /2m, ip = ip + ^rSgn(p), and /i is the 
chemical potential. In order to obtain x i n the real frequency, analytic continuation of iu — ► cj + i<5 is performed after 
the frequency summation. The frequency summation can be done in complex plane using contour integral, where 
summation is replaced by inegrals over branch cut lines (q| . Then, Ximp can be written as 



Xim P (.q,w)=-i del [n F (e + u>) - n F (ej\ F(e - — ,e + w+ — ) 
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where we introduce for covenience an expression 
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and 
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Here, n F (e) = 1/[1 + exp(/3e)] is the Fermi distribution function. 

Now we can express K and (TiT) -1 with this susceptibility, through taking its real and imaginary parts. (XiT) -1 
can be written as 
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where X2 and Ximp2 denote the imaginary parts of \ and Ximp: respectively. Ximpi and Ximp2 in Eq. (g) can be written 



as 
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where B = /3e F , r\ = (2re F ) , M = /J,/e F , t = q/k F , x = s/e F , and 



T , , . . . i , (±\fx+M±iii+y/x+M + iri + 1 

J(x±iri,x+rn) = — In = = 

" It \±^/x+M±ir)+^x+M + ir)-t 



(10) 
(11) 

(12) 



In the same way, K is written as 
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where Ximpi(0, 0) is given from Eq. ( [io| ) by 
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r\J in the denominater in Eqs. ( |10[ ) and (|ll|), and inoutside the sqare root in Eq. ( J14| ) are from the impurity vertex 
correction. These vertex parts from Eqs. fllO| ) and ( |l4| ) cause the crosses between the curves of K. with T and [/. 



III. RESULTS 



Let us first consider the free electron gas model for a reference before we present the calculations for ferromagnetically 
correlated dirty metals. For a 3D free electron gas model, it is simple to calculate the static susceptibility x(q, w — > 0) 
as a function of T M. From x, K and (TiT) -1 are given as 
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where the chemical potential /j is given by 
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We, therefore, have 
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where i^To = 2^%^ and (TiT) x = 2tt/i 2 3 N 2 denote, respectively, their zero temperature values at U — 1/r = 0. From 
these elementary considerations, we see that 
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Introducing impurity scatterings and electron-electron interaction does not alter in an essential way the K, \jT\T 
and K, vs. T behavior of Eqs. (p"7J), (pf), and (Il9|), but change their T dependence by altering the three coefficients: 
(a) the value at T = 0, (b) the low T coefficient proportional to T 2 (Lorenzian width), and (c) the high T coefficient, 
as shown in Figs. 2 and 3. K/Kq as a function of T for a set of r\ = l/(2re F ) and NU, where NU cr = 1 is the critical 
value at which the ferromagnetic instability occurs at T = and ?y = 0, are shown in Fig. 2(a). Both NU and r\ enhance 
K in agreement with Shastry and Abrahams. The low temperature coefficients of K/Kq vs. T can be calculated 
analytically as functions of NU and r\: (a) The T — peak is given by a/(l — aNU), and (b) the low T Lorentzian 
width by 2a/I - aNU/Vb, where a = 1 + f 77 - fr? 2 and 6 = f|[l - (1 + f )?? + ±(3 + n + ^)t] 2 }. The numerical 



calulations agree perfectly well with the analytically obtained asymptotic behavior. In Fig. 2(b), (TiT) _1 /(TiT) _1 as 
a function of T are shown for a set of 77 and NU. Their general behavior is almost the same with that of K/Kq. In Ref. 
H and p), (TiT)^ 1 measured as a function of T is almost constant above T c in the fullerene superconductors. It is 
difficult, however, to compare the present results with experimental observations because the lattice constant changes 
from the thermal expansion. The lattice constant expansion affects the DOS and, consequently, K and 1/T\T. 

In Fig. 3(a), we show K. as a function of T for a set of ij and NU. One can see that (a) the T — value is suppressed 
by NU, but, by i], it is increased for small U and decreased for large U, (K = 1 for T = U = r) = 0.) (b) the low 
T coefficient is nearly independent of NU and 77, and (c) the high T coefficient proportional to T 3 / 2 is reduced as rj 
becomes larger. Above changes cause the K. vs. T curves cross each other as can be seen in Fig. 4. These crosses do 
not appear if we neglect the impurity vertex correction in Ximpi, that is, if we neglect the 1/r term in the denominator 
of Ximpi, and their origin is this impurity vertex correction, as will be discussed in more detail below. Now, let us 
consider the T = limit case, which can directly be compared with Shastry and Abrahams. In Fig. 3(b), we show /C 
as a function of NU for r\ = 1CP 4 , 0.05, and 0.1 at T = 0. The three K, vs. NU curves cross each other as NU is 
increased, as in Fig. (a) discussed above. These crosses did not occur in Shastry and Abrahams, who neglected the 



impurity vertex correction in Ximgi- At T = 0, 
be expanded by integrating Eq. (|6|) as, 



and for small q and u), the impurity enhanced susceptibility, Ximp, can 
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where the diffusion constant is given by D = ru F /3, and Xi mp (0; 0) = Xi mp (0, 0) = Xim P i(0, 0) = 1 + \ 77 + |?7 2 , up to 
the order of ?7 2 . The 7r?7/4 in Xi mp comes from the impurity vertex correction in the u> = q = limit, and ?7 2 /8 from 



the self-energy correction due to impurity scatterings in the renormalized Green's funtions. Shastry and Abrahams 
neglected the impurity vertex corrections in calculating Ximp(<7, w = 0), and found that x im (0,0) sa 1 + f] 2 /8 in 
agreement with the present result. These two terms enhance Ximpi(0,0) compared with the free electron value of 1. 
Because Ximpi(<Z, 0) is multiplied by NU in the denomenator when calculating K or 1/(T\T) as in Eqs. (||) or (|l3|), the 
effects of U are stronger for larger rj. Consequently, K decreases faster for larger r\ as NU is increased. This, together 
with the fact that K. is larger for larger rj at U = 0, causes the crosses among the curves of K. vs. NU as shown in 
Fig. 3(b) Since the 7rr;/4 term in Ximpi(0, 0) is dominant contribution from the impurity scatterings compared with 
r] 2 /8, these crosses do not occur without the impurity vertex correction. It should be interesting if this prediction can 
be confirmed by experiments. 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, we studied the Korringa ratio K. = l/(TiTK 2 ), Knight shift K, and nuclear spin-lattice relaxation 
time T\ of ferromagnetically correlated impure metals at finite temperature. They were calculated from the spin 
susceptibility including the self-energy correction due to the impurity scatterings, and the vertex correction due 
to the 6- function electron-electron interaction and impurity scatterings. The spin susceptibility was formulated with 
Matsubara Green's function to consider the finite temperature effects and analytically continued to the real frequency. 
The K and (TiT) -1 decrease with T, but K. increases as T is raised because K 2 decreases more rapidly than (TiT) -1 
with T. Also, K and (TiT) -1 increase as U or rj is increased, and K. decreases as U is increased. 

The conventional interpretation of the Korringa ratio, as mentioned in Introduction, is that when K of a material 
is larger (smaller) than unity, it implies that the material is antiferromagnetically (ferromagnetically) correlated. The 
Shastry and Abrahams contribution is that K. is increased by the non-magnetic impurity scatterings, so that K. of 
ferromagnetically correlated dirty metals, where K. is suppressed by the ferromagnetic correlation but increased by 
the impurity scatterings, can be close to 1, just like the free electron case. This means that one should be careful in 
extracting the strength of ferromagnetic correlation from the Korringa ratio. 

Our observation in the present work is that the Korringa ratio is increased by the non-magnetic impurity scatterings 
when both U and T are small, but K is decreased by the impurity scatterings when U is comparable with U cr or T 
is comparable with ej?. This implies that K is enhanced by the impurity scatterings only for weakly ferromagnetic 
dirty metals at low T, but supressed when U or T is large. This causes the crosses among the K vs. NU curves for 
the different impurity scattering rates 1/r in the present calculations. The crosses also occur in K. vs. T at a given 
U among diffcrncnt values of 1/r. Those crosses are caused by the impurity vertex correction in Ximpi(0,0), which 
was neglected in the previous calculations. The present exact numerical calculations within the finite temperature 
extension of the Fulde and Luther formulation, shown in Fig. 3, confirm that above behavior is due to the impurity 
vertex corrections. This observation requires that the interpretation according to Shastry and Abrahams should be 
modified, and one should do more careful analysis in estimating the strength of ferromagnetic correlation of impure 
materials from the Korringa ratio. 
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FIG. 1. (a) is the Feynman diagram of the spin susceptibility, and (b) is for the renormalized vertex F including the 
impurity and Coulomb interactions. The solid line indicates the electron Green's function, the dashed line the electron-electron 
interaction, and the dashed line with t a cross stands for the impurity scatterings, ip and iu> are, respectively, the fermion and 
boson Matsubara frequencies, and k, q are correspinding wave vectors. 
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FIG. 2. (a) is the normalized Kight shift K/K vs. the reduced temperature, T/e F , and (b) is (TiT) _1 /( T i T )o X vs. T/e F . 
Kq and (TiT)o - are, respectively, the zero temperature values of K and (TiT) -1 at T — U = r\ — limit. NU = 0, 0.4, 0.8 
from bottom to top and ef = 0.25 eV. K/Ko and (TiT) _1 /(TiT)(P decrease more rapidly with T when NU or r\ is large. 



K 




(b) 


— o— 77 = 1 0" 4 




— . — n = 0.05 




77 = 0.1 




x\ 








XSx A 




"^ 


s, 








\ 


- 


V 



0.0 



0.2 



0.4 



0.6 



1.0 



T/e c 



NU 



FIG. 3. (a) is the Korringa ratio K, vs. T/e F , and (b) is K, vs. NU at T = limit. £=1 in the T = (7 = ry = limit. 
ef = 0.25 eV, and NU — 0, 0.4, 0.8 from bottom to top for Fig. (a). The crosses among the curves for different 77 occur at lower 
T as NU is increased, implying that in K. vs. NU crosses occur at lower NU when T is increased. Fig. (b) can be compared 
with Ref. 3, where there are no crosses among different r\ values. 



